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Abstract 

Q-| The i? 7^ decay at large energies of the photon receives a numerically impor- 

tant soft-overlap contribution which is formally of the next-to-leading order in the 
expansion in the inverse photon energy. We point out that this contribution can be 
calculated within the framework of heavy-quark expansion and soft-collinear effective 
theory, making use of dispersion relations and quark-hadron duality. The soft-overlap 
contribution is obtained in a full analogy with the similar contribution to the 7*7 — > n 
transition form factor. This result strengthens the case for using the B — > ^ivg decay 
to constrain the 5-meson distribution amplitude and determine its most important 



parameter, the inverse moment A#. 



1. The decay B — > j£u£ at large photon energy is one of the simplest hadronic 
processes that can be studied using QCD factorization [TJ |2j [31 H] . This method involves the 
B-meson light-cone distribution amplitude (DA) [SJ El H El El EE] as the main nonperturbative 
input at the leading order in the heavy quark expansion. The B — > 7^ decay U is therefore 
best suited for determining the parameters of the 5-meson DA and presents a close analog to 
the photon-pion transition 7*7 — > it with one highly-virtual {Q 2 ) and one real photon, which is 
used for the determination of the pion DA. The theoretical challenges are also similar, and in 
both cases are mainly due to the necessity to have a quantitative control over the corrections 
that are subleading in the expansion parameter, l/m b ~ 1/(2E^) or l/Q 2 , and, in general, 
are not factorizable. The purpose of this letter is to emphasize this connection and suggest 
a method to calculate the soft contribution to B — > 7^, which closely follows the technique 
used in 7*7 — > n. 

A recent summary of the theory status of B — > "jlvg can be found in [TTJ. The decay 
amplitude 

A(B~ j£u e ) = ^^(£un\£Y(l - 75) - l5)b\B~) (1) 



Ho distinguish this important decay channel we suggest to call it photoleptonic -B-decay. 
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can be written in terms of the two form factors, Fy and Fa, defined through the Lorentz 
decomposition of the hadronic tensor 

T^( P ,q) = -tjd'xe^(0\T{ 3 ; m (x)miAl-l 5 )m}\B-(p + q)} 

= e^ Tp p T v p F v +i[- g^ipv) + v^] F A + ... . (2) 

Here p and q are the photon and lepton-pair momenta, respectively, so that (p + q) = vtlbv 
is the 5-meson momentum in terms of its velocity. In the above, j£ m = J2 g ^QlnQ is the 
electromagnetic current and the ellipses stand for the terms ~ p^ and for the contact term. 
The origin of the latter is explained in [TT] (see also p2]). The form factors can be written as 
functions of the lepton-pair invariant mass squared q 2 , or, equivalently, of the photon energy 
Ery in the 5-meson rest frame. The differential decay width is given by 

dT _ ^Gl\V ub \ 2 f 1 _^f| fy | 2+ \f a A , F A = Fa + —Jr~ , (3) 



dEj 67r 2 7 y m B J ^ ' E 1 

where the contact term is included in the axial form factor, and the lepton mass is neglected. 
For large photon energies the form factors can be calculated [11] as 



e u fBm B 



ebfBm B , e u f B m B 



Y2 



p.pN e u f B m B \cftr^ e b f B m B e u f B m B ] 

The first term in both expressions represents the leading contribution in the heavy-quark 
expansion (HQE) that corresponds to the photon emission from the light spectator quark in 
B meson. In the above, f B is the decay constant of B meson, and the quantity As is the first 
inverse moment of the 5-meson DA: 

oo 

— (t>+(u,fj), (5) 







where the variable u is the energy of the light-quark in the 5-meson (in its rest frame). The 
factor R(E^,/j,) in (TjJ takes into account gluon radiative corrections (see [TT] for details) and 
equals one at the tree level. The scale-dependence of + (w,/i) was calculated to the leading- 
logarithmic accuracy in [T3] . 

The terms in square brackets in flj]) are the l/m& and 1/(2E 7 ) power corrections to the 
leading-order expression. They are written in this particular form to emphasize that some of 
these corrections are "symmetry-preserving" , i.e. are the same for both form factors Fy and 
Fa, and some corrections are "symmetry-breaking", i.e. they are different. The symmetry- 
preserving corrections parametrized by the function £(i£ 7 ) present the main difficulty They 
are unknown, apart that from the power counting one expects £(i? 7 ) ~ I/-E7 with respect to 
the leading-order term. In the analysis of [11 J this function was parametrized as 

£(£ 7 ) = cJ-, (6) 
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that is, tacitly assuming al/ra^ suppression with respect to the leading-order term in (jlj), and 
the coefficient c was varied between —1 and +1. 

In this letter we present an approach which allows one to estimate £(i? 7 ) to, potentially, 
20% accuracy. To be precise, we will be calculating the soft overlap contribution to £(E^) 
which is not directly accessible in QCD factorization. Generally speaking, there exist also 
factorizable symmetry-preserving contributions which can treated in a systematic way within 
the HQE and added. 

2. The main idea which parallels the technique originally suggested in [TJ] for 7*7 — > tt 
transition form factor is to consider the hadronic tensor (T5j) at negative p 2 < 0, m 2 B ^> \p 2 \ ^> 
Aq CD , or in other words, an unphysical decay B — > j*£v£ involving a (transversely polarized) 
spacelike photon. The corresponding form factors, now functions of two variables q 2 and p 2 , 
can be calculated in QCD using HQE and operator-product expansion (OPE) to (at least in 
principle) arbitrary power accuracy in l/m^, 1/E^, 1/p 2 . The analytic continuation to the real 
photon limit p 2 = can be made using dispersion relation. In this way the explicit evaluation 
of soft nonfactorizable contributions is effectively replaced by a certain ansatz of the hadronic 
spectral density in the p 2 -channel. 

The starting observation (cf. [H]) is that Fv,a{q 2 ?P 2 )> the generalized form factors of B — > 
7*£^, at fixed q 2 satisfy an unsubtracted dispersion relation in the variable p 2 . Separating the 
contribution of the lowest-lying vector mesons p, u>, one can write (F B -> 7 * = Fy or Fa) 



(9 V) = />wf) + 1 lis t (7) 



m 2 — p 2 Ti J s — p 2 

so 



where s is a certain effective threshold. Here, we combined the p and u contributions in 
one resonance term, assuming m p ~ m u and adopting the zero-width approximation. In the 
above, f p is the usual decay constant of the vector meson and F B ^ p (q 2 ) is a generic B — » p{uj) 
transition form factor. Note that since there are no massless states, the real photon limit is 
recovered by the simple substitution p 2 — >■ in (JTj). 

On the other hand, the same form factors can be calculated using HQE. The result satisfies 
a similar dispersion relation 



1 7 ImF HQE (a 2 s) 



TT J S — £T 



where the limit p 2 — > cannot be taken directly: e.g., singular terms in 1/p 2 appear (cf. |15j). 
signaling that the HQE cannot be applied to the real photon case p 2 = beyond the leading 
power in I /nib. 

The main assumption of the method is that the physical spectral density above the threshold 
so coincides with the QCD spectral density, as given by the HQE: 

lmF B ^ r (q 2 ,s) ~ImF^(g 2 , S ) for s > s . (9) 

This is the usual approximation of local quark-hadron duality. In reality we employ a weaker 
form of duality, assuming that ([9]) holds upon averaging with a smooth weight function over a 
sufficiently broad interval of the variable s. 
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We expect that the HQE reproduces the "true" form factors Fb_>. 7 * (q 2 , p 2 ) for asymptot- 
ically large values of — p 2 . Equating the two representations (|7j) and (JHJ) at p 2 — > — oo and 
subtracting the contributions of s > so from both sides one obtains 

SO 

f p F B ^{q 2 ) = lj dslmF^(q 2 ,s) . (10) 
o 

This relation explains why sq is usually referred to as the interval of duality in the vector- 
meson channel. The perturbatively obtained spectral density in HQE ImF B _*j^.(q 2 , s) is a 
smooth function of s and does not vanish at s — > 0. It is very different from the hadronic 
spectral density lmF B ^*£ Ue (q 2 , s) ~ S(s — m 2 ). However, the integrals over both spectral 
densities over a certain region of s coincide; in this sense the QCD description of correlation 
functions in terms of quark and gluons is dual to the one in terms of hadronic states. 

In practical applications of this method one uses an additional device, borrowed from the 
method of QCD sum rules [TB], which allows one to reduce the sensitivity to the duality 
assumption in and simultaneously suppress the contributions of higher orders in the OPE. 
This is done going over to the Borel representation l/(s — p 2 ) — > exp[— s/M 2 ], the net effect 
being the appearance of an additional weight factor under the integral 

SO 

f P F B ^ P {q 2 ) = lj dse-^/ M2 lmF^(q 2 ,s) . (11) 
o 

Varying the Borel parameter within a certain window, usually M 2 = 1 — 2 GeV 2 , one can test 
the sensitivity of the results to the particular model of the spectral density H. 

With this refinement, substituting ffTTl) in ([7]) and using Q one obtains for p 2 — >■ (cf. [TH 

US]) 




This expression contains two nonperturbative parameters — the vector meson mass m p and 
effective threshold so — 1.5 GeV 2 — as compared to the "pure" QCD calculation based on the 
HQE expansion. The reward is that the HQE can be done to an arbitrary accuracy in the 
powers of l/m^. The nonfactorizable contributions that are beyond the accuracy of the HQE 
in the usual treatment are taken into account effectively, via the nonperturbative modification 
of the spectral density. 

3. As an illustration, consider the expression corresponding to the leading-order in 
HQE diagram of the virtual photon emission from the light spectator-antiquark in B meson. 
Calculating this diagram in terms of 5-meson DA at p 2 < 0, we obtain, after replacing q 2 = 
m 2 B — 2m B E^ +p 2 , the following expression for the form factors defined in ([2]): 

oo 

F$\q 2 ,p 2 ) = F B °\ r (E„p 2 ) = e u f B m B [ du ^ U ^\ , (13) 



2 We note in passing that the relation where the r.h.s. is calculated in terms of HQE and B meson 
DA's, represents a light-sone sum rule for B — ^ p form factor of the particular type considered in [17| 118] . 
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where we neglected the corrections ~ u/E 1 . In this approximation, the axial form factor 



x A (q 2 ,p 2 ) coincides with F v '(q 2 ,p 2 ). The energy integral in (|T3|) can easily be converted to 
the form of a dispersion relation by the change of variables s = 2E^u. Following the procedure 
described above and changing the integration variable back to u 
p 2 0, 



s/(2E y ), we obtain, at 



? (o) 



(£ 7 ) 



IE™ 



so/(2B 7 ) 



J rni J u 



(f>+(u,(j,) 



(14) 



s /(2£ 7 ) 



Completing the second integral to run from zero to infinity and subtracting the correction from 
the first term, we can write this expression as 



7 (0) 
B^-y 



(E, 



e u fBm B , e u f B m B 



s /{2E 1 ) 



+ 



2E 1 \ B (ii) 2E, 



du 



2 -^7 c ~(2E^u-m 2 r )/M 2 _ }_ 



mi 



U 



(15) 



where the first term is nothing but the HQE expression for the form factor to leading order in a s 
and leading power accuracy in (j2J). The second term can be identified with the soft correction 
^(E^) as it appears in the same equations. Note that the modification of the standard HQE 
expression only concerns the region u < sq/{2E 1 ) ~ So/rrib, hence, this is a soft spectator-quark 
contribution. If <f) + (u,fi) ~ u forw^OjgiHEllHJCEn], C (0) (^7) defined by $H]) corresponds 
to a power correction of the order of Sq/{2E 1 ) 2 for E 1 ~ — > oo, in agreement with usual 
power counting. 

We define a rescaled soft contribution £(E~) such that 



(E 7 ) 



' e u f B m B 



+ . . . 



(16) 



where the expression in the parenthesis is the leading-order result, and we anticipate that the 
function £(E^) depends on E 1 only weakly. The ellipses stand for radiative and "hard" power 
corrections, cf. (jl]) %. In the adopted approximation 



so/(2B 7 ) 



^\E 1 ) = 2E 1 \ B J 



du 



2E, 



U 



;i7) 



Note that the definition of the S-meson DA in the soft-collinear effective theory (SCET) 
involves a collinear light spectator-antiquark field. If the separation between collinear and soft 
regions were done with an explicit cutoff, u > /j,g C /E 7 , the leading-order result for the form 
factors would read 

oo 

e u fBm B f du 



w SCET ( -p \ 

i's^ {E~ f/ 



u 



3 E.g., the corrections given by the ~ e u terms in square brackets in ((4|) are readily obtained if one retains 
the 0(u!/E-y) terms in the integral (fT3|) . 
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Figure 1: The soft contribution £(22 7 ) to the form factors in B — > j£u£ (left panel) for the first (solid 
curves) and second (dashed curves) model of the 5-meson DA in (I19p , compared to the soft contribution 
£{Q 2 ) to the 7*7 — > tt form factor (right panel) for the asymptotic pion DA (solid) and a realistic model 
with 02,4 7^ (dashed). The shaded areas correspond to the variation of the Borel parameter in the range 
M 2 = 1.0-1.5 GeV 2 . 



It is easy to see that the soft correction defined by f fT4|) . f lT5|) effectively cuts off the small 
energy region in a way similar in spirit to ( IT8|) . with the interval of duality playing the role 
of the hard cutoff /isc- We therefore expect that the soft nonperturbative correction is always 
negative relative to the leading-order result because its role is, conceptually, to create a mass 
gap in the vector-meson mass spectrum. 

4. For numerical estimates we will use two models of the S-meson DA: 



0+M 



u 

A 2 "' 



-ui/\b 



ttX b UJ 2 + 1 



1 



2{<T B - 1) 



u 2 + l 



In Co 



7l z 



(19) 



where Cj = oj/1 GeV and gb is the first logarithmic moment ctb\ b =—f — 0+(w,/i) In-. We 

o M 

take 300 MeV < X B < 600 MeV and a B = 1.5 [HI E] (see also [TO]) as typical values of the 
parameters. 

The function £(£7) calculated for the two models of the B- meson DA assuming the standard 
choice of the continuum threshold s = 1.5 GeV 2 and M 2 = 1.0 — 1.5 GeV 2 is plotted vs. 2E 7 
in Fig. [U the left panel. For this plot we have chosen A# = 500 MeV. Note that the curves are 
essentially flat for 2E 1 > 1 GeV which means that the soft contribution is well described by a 
power-suppressed correction ~ l/(2E y ) as compared to the leading-order result. 

For comparison, we show in the same figure (right panel) the soft correction to the 7*7 — > 71 
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transition form factor in the same approximation: 



V2fn f f dx 




X(, 



k(x) + dx 



Q 2 e (xm 2 p ~xQ 2 )/(xM 2 ) _ }_ 



xm? p 



x 



—(f) w {x) 

X 



£7*7— >n{Q ) 

Q 2 



(20) 



where xo = so/(so + Q 2 )- The results are shown for two sample models of the pion DA: 1) 
the asymptotic DA (f> % (x) = 6x(l — x) and 2) a realistic DA with nonvanishing Gegenbauer 
moments a,2(lGeV) = 0.17 and a^lGeV) = 0.06. These are typical numbers that are used in 
the light-cone sum rule analysis of the pion electromagnetic form factor [19] and weak B — > wife 
decays [20], see also a discussion in [T5] . 

Note that the size of the soft correction to the 7*7 — » tt form factor for the asymptotic 
pion DA is very similar to the soft correction to B — > 7^ for the existing models (and for 
Xb = 500 MeV). For the realistic pion DA the correction is larger. The difference is partially 
due to the larger value of the inverse moment / (dx/x)<p n (x) = 3(1 + a 2 + a 4 + . . .) playing the 
same role as A^ 1 for 5-meson, but also due to the functional form. If the pion DA has some 
enhancement close to the end points, as was suspected in particular in view of the BABAR 
data [21] but also not excluded by Belle (22], then the soft correction is larger and the onset 
of the asymptotic regime (where it is a 1/Q A correction) occurs much later, see [T5| 123]. 

Varying the inverse moment of the 5-meson DA in the interval 300 MeV < Xb < 600 MeV 
we obtain the following values of the rescaled soft factor £, defined in (flE]) . averaged over the 
photon energy interval 2 GeV < 2E^ < tub'- 



X b [MeV] 


300 


400 


500 


600 


T^T 

£ Model I [GeV] 


- (0.5018:3) 


- (0.3618:2) 


- (o.27l8:8 9 7 ) 


- (0.2218:8^) 


(.Model II [GeV] 






- (0.2318:8^) 





The quoted uncertainties include variations of the Borel mass and the photon energy depen- 
dence. 

Our result for £ (Model I) translates to the value of the coefficient c in the notation of ([6]): 



m B 
2EL 



a/ 



m B 
2E^ 



x < 



-0.09\ 
-0.27J > 

(o.6oi8:l°), 

-0.09\ 
-0.12J 1 

(0.2418:8^), 



'l.ir 



-(0.36] 



X B 
X b 
X b 
X b 



300 MeV 
400 MeV 
500 MeV 
600 MeV 



(21) 



where the uncertainties have the same origin as in the above Table. We emphasize that the soft 
contribution £(E 7 ) obtained here has a suppression factor 1/(2E 7 ) with respect to the leading 
term, as compared to the 1/uib suppression assumed in [UJ. Note that for large A^ 1 , i.e. for 
the 5-meson DA that is enhanced in the soft region, there is a strong cancellation between 
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Figure 2: The vector and axial form factors for B — > jiv^ calculated from (TJJ at A# = 500 MeV. 
The lowest curves show the soft-overlap part £(E^) of the form factors. The shaded areas indicate the 
uncertainties estimated by varying the Borel parameter within M 2 = 1.0 — 1.5 GeV 2 and switching from 
the model I to model II of the S-meson DA. 

the leading term and the soft contribution: In a hypothetical limit A# — > both terms diverge 
but the sum of them remains finite. 

5. We can employ our estimate of the soft-overlap contribution to calculate the form factors 
Fy and Fa at large photon energies. To this end we use the expressions in (j3J) and the result 
for R(Ery, /i) given in [TT] which includes resummation of radiative corrections to the next-to- 
leading order logarithmic accuracy. Following [11 J we adopt the "soft-collinear" scale /x = 1.5 
GeV, the heavy-quark mass scale = 4.8 GeV and /# = 195 MeV. In the context of this 
study we are interested mostly in the uncertainty due to the soft contribution; the variation 
of the i?-meson decay constant within the intervals of the current lattice QCD estimates, will 
yield an additional (correlated) theoretical uncertainty in both form factors about ±10%. The 
results for the vector and axial form factors — the latter including the contact term — and, 
separately, for the soft contribution defined in our approximation as f (£ 7 ) = f^l^ £ (0) (ff 7 ), 
are shown for the choice Xb = 500 MeV in Fig.PJJ From fl2T|) it is clear that our estimate of 
the soft form factor has a smaller error than the interval — 1 < c < 1 taken in |llj . 

Using these form factors we calculate the partial branching fraction BR(B — > 'jive), inte- 
grating ([3]) over the photon energies E min < E 7 < m B /2 . The result is shown in Fig. [3] as a 
function of A# for two different choices of the photon energy cut, E min = 1.0 (1.7) GeV. For 
this plot, for definiteness, we take \V u b\ = 3.5 x 10~ 3 , in the ballpark of current determinations 
from exclusive semileptonic B decays. Our main message is that the uncertainty due to the 
soft overlap contribution is sufficiently small. 

The existing experimental data are not yet conclusive. The upper limit on the partial 
branching fraction BR(B —>• jivt) with certain kinematical cuts, including the cut on E 7 , was 
published by the BABAR collaboration in [24J. This limit is weaker than their previous result 
quoted in [25]. As explained in [11], the published limit [23] is not yet sufficient to constrain 
the inverse moment A#. 

Finally, we suggest to consider the ratio of the photoleptonic and leptonic charged I?-meson 
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Figure 3: The partial branching fraction BR{B — > "f£vi)E 1 >E min for E m i n = 1.0 GeV (upper, solid) and 
Emin = 1.7 GeV (lower, dashed). The uncertainties indicated by shaded areas are of the same origin as in 

Fig. m 



decay widths: 

BR(B ->■ 7M^) E „ >BmiK 



BR(B -»■ rz/ r 



m B /2 



lM "" W*fl — 



37rm 2 (l — m 2 /m 2 B ) 2 J 7 7 ^ 



|F y (E 7 )| 2 + |F A (E, 



,(22) 



where we neglect the muon mass. Both observables are accessible in the 5-factory exper- 
iments and their ratio does not depend on V U b and on the 5-meson decay constant be- 
cause Jb enters also the normalization of the form factors. We predict: R lilv lTv{E m %n = 
1.7 GeV) = 0.0103,0.0058,0.0037,0.0025 at A B = 300,400,500,600 MeV, respectively (for 
model I, M 2 =1.0 GeV 2 ). Note that the recent measurement [26] of B — > tu t yields /b|1^&| = 
[7.4 ± O.S(stat) ±0.5(syst)} x 10~ 4 GeV, consistent with the lattice-QCD value of and with 
the value of V u b which we have used above for the estimate of the partial photoleptonic B 
decay width. 

6. To conclude, we have described a method to calculate the soft contribution to the 
decay B — > ^Iv^ which is formally subleading in powers of the photon energy E 1 and is 
not directly accessible in QCD factorization. This method has originally been developed for 
another process, the 7*7 — > n transition form factor, in which case the QCD calculation is 
under a better control because the moments of pion DA can be calculated in lattice QCD. The 
successful description of the experimental data on the 7*7 — > n form factor in this framework 
in the region of momentum transfers Q 2 = 2 — 6 GeV 2 (see e.g. [13] and references therein) 
allows one to hope that the same technique will yield sufficiently accurate predictions for the 
decay B — > j£i/£ as well. 

The calculation presented in this letter serves the purpose of demonstration mainly. It can 
and should be extended in several directions, by replacing the leading-order expression for the 
spectral density ImF^^, in ( TT2"|) by the complete HQE expression to the 0(l/(2E / y), l/m B ) 
accuracy. First of all, one has to include radiative corrections which give rise to the R^E^,^) 
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factor in This requires a calculation of the terms 0(a s ) in the coefficient function of the 
S-meson DA for nonzero photon virtualities, which is straightforward. Such corrections can 
then be resummed using SCET techniques, although we expect that the numerical impact of 
the resummation will be minimal. 

Second, one has to include contributions of higher-twist two-particle and also three-particle 
S-meson DAs. These contributions contain logarithmic end-point singularities if calculated 
directly, but give rise to finite contributions to the dispersion integral ( fl2|) with the continuum 
threshold So providing an effective IR cutoff. In other words, the full contribution of twist- 
four DAs to the form factor is finite, but an attempt to rewrite the answer as a sum of the 
"pure" HQE expression plus a correction, as in ( TT51) would result in divergent expressions. We 
believe that the corresponding calculation would be very interesting because the normalization 
constants in the higher-twist 5-meson DAs are related to higher moments of the leading-twist 
DA, at least within certain schemes used to subtract the corresponding high-energy behavior, 
see [5J mHOl [27]. 

Third, there are terms related to photon emission at large distances that involve a photon 
DA. These contributions were estimated for the 7*7 — > it transition form factor in [TB] in which 
case they correspond to contributions of twist-six operators in the operator product expansion. 
This example is instructive as it shows that for soft corrections the correspondence between 
power suppression and twist counting is lost: Contributions of all twists to the operator product 
expansion of the electromagnetic form factors produce power corrections which are suppressed 
by one power of Q 2 with respect to the leading twist result. We expect that the situation with 
the HQE for the decay B — > ^tvi is similar. 

Last but not least, let us mention that B — > jiisg decay amplitudes at large photon en- 
ergies have also been calculated using light-cone sum rules with photon DA's and S-meson 
interpolating current [2"8"| I2"9~t 150"] . It would be interesting to clarify the interconnection of this 
approach with the methods employing the 5-meson DA's in order to gain a more complete 
picture of the underlying quark-gluon dynamics. 
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